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METHOD AND APPARATUS FOR NEUTRON MICROSCOPY WITH 

STOICHIOMETRIC IMAGING 

Related Applications 

The present application claims priority to U.S. Provisional Patent Application No. 
60/213,373, entitled "Method And Apparatus For Neutron Microscopy With Stoichiometric 
Imaging" filed June 23, 2000. The present application also claims priority to U.S. Patent 
Application No. 09/788,736, entitled "Method And Apparatus For Detecting, Locating, And 
Analyzing Chemical Compounds Using Subatomic Particle Activation" filed February 20, 2001, 
which is a continuation from U.S. Patent Application No. 09/265,043, entitled "Method And 
Apparatus For Detecting, Locating, And Analyzing Chemical Compounds Using Subatomic 
Particle Activation" filed March 9, 1999, which is a continuation-in-part from U.S. Patent 
Application No. 09/252,359, entitled "Method And Apparatus For Detecting, Locating, And 
Analyzing Chemical Compounds Using Subatomic Particle Activation" filed February 17, 1999, 
which claims priority to U.S. Provisional Patent Application No. 60/075,037, entitled "Method And 
Apparatus For Detecting, Locating, And Analyzing Chemical Compounds Using Subatomic 
Particle Activation" filed February 18, 1998. 

Background of the Invention 

1. Field of the In vention 

The present invention relates to the field of non-invasive detection, imaging and 
stoichiometric analysis of inaccessible minute or microscopic quantities of chemical compounds 
and elements, and of mechanical structures and flaws inside materials, using subatomic particle 
activation. Here, stoichiometric means deciphering of the empirical chemical formulas of 
substances. 

2. Description of Related Technology 

The manifold societal problems associated with unexploded antipersonal ("AP")/anti- 
vehicular land mines and chemical weapons are well known and documented. Such problems 
include, inter alia, the inadvertent detonation of such devices by an unsuspecting civilian population 



often times many years after the cease of hostilities. As of the late twentieth century, vast portions 
of the surface of the earth are infested with such devices and therefore rendered largely unusable. 

In addition to land mines, significant stockpiles of unexploded chemical weapons in the 
form of artillery shells, rockets, grenades, and other warheads exist. These weapons contain a 
variety of highly destructive and potentially lethal compounds such as Sarin, and often bear no 
markings or means of identification of their contents thereby making disposal highly inefficient and 
dangerous. 

Furthermore, large amounts of chemical contraband (e.g., illicit drugs such as cocaine, 
heroin, marijuana, and PCP, or alcohol) are produced and distributed throughout the world on a 
daily basis. These substances result in a host of deleterious effects on society in general including 
increased health care and rehabilitation costs as well as constant monitoring, surveillance, and 
intervention by law enforcement agencies. 

Existing techniques for detecting, physically locating, and analyzing the aforementioned 
chemical compounds are marginally effective at best. A variety of different techniques such as X- 
ray analysis, magnetic resonance imaging (MRI), chemical "sniffers", and visual inspection have 
been employed to date, yet all suffer from one significant disability or another, thereby greatly 
reducing their efficacy. For example, X-ray techniques can only provide information about an 
object's shape or location, and are not useful in large area searches (such as for land mines buried in 
the field or searches of large containerized cargo). Furthermore, such techniques require the 
subsequent use of intrusive means to determine if the identified substance is dangerous or not, thus 
resulting in a very high proportion of "false alarms." Chemical sniffers are effective under certain 
limited circumstances, but can be easily defeated through proper sealing of the chemical compound 
in a non-permeable container, and are also impractical for use in many applications. Recently, more 
promising methods of detection and analysis using nuclear radiation (including so-called "fast 
neutron activation" or FNA techniques, such as described in U.S. Patent No. 5,098,640, "Apparatus 
and Method for Detecting Contraband Using Fast Neutron Activation") have been developed, yet 
these methods still suffer from a number of problems of their own, including poor spatial and 
gamma ray spectral resolution, great size, weight, and complexity. One significant problem related 
to these systems concerns the use of prior art subatomic particle coincidence circuits, which operate 
on the principle of detector-to-detector (or counter-to-counter) coincidence. This approach 
necessitates an analysis of the entire gamma spectrum generated by the counter, thereby requiring a 



-2- 



tremendous signal processing capability to analyze even a modest number of detection events per 
unit time. The net result is very long irradiation/counting times (and correspondingly lower incident 
particle flux), as well as reduced chemical identification accuracy and confidence. Additionally, 
poor energy resolution of the scintillation detectors used in these systems has hampered the 
identification of specific spectral artifacts. These disabilities are discussed in greater detail in the 
following paragraphs. 

Chemical Detection and Identification Using FNA-Induced Gamma Rays 

There are two primary requirements for the quantitative chemical detection of explosives 
and other contraband substances using gamma rays: (i) the ability to resolve gamma ray energy 
precisely (approximately 0.5% or better resolution required); and (ii) the ability to temporally 
resolve gamma ray events (approximately 3 nanosecond and better resolution required). Prior art 
FNA devices have characteristically opted for good time resolution (typically 1-2 ns), at the 
expense of energy resolution. For example, a typical prior art Sodium Iodide scintillation 
detector has an energy resolution on the order of 10%. Because of this comparatively poor 
energy resolution, prior art contraband detection devices based on gamma ray spectrum analysis 
are not able to determine quantitative elemental composition of the interrogated object. Rather, 
they can detect only the presence (or dominance) of certain elements that are potentially 
indicative of contraband, thus providing an 'alert' signal. This alert signal in turn requires 
intrusive inspection of the interrogated object, and results in a false alarm rate which, while 
reduced from that of X ray contraband detectors, is still quite high (over 90% by some 
estimates). In contrast to scintillation detectors, High Purity Germanium Detectors (HPGD) have 
an energy resolution on the order of 0.1-0.3%; however, their dead time, determined by their 
charge collection time, is roughly 200 ns. Such dead time was considered incompatible with the 
aforementioned requirement of high temporal resolution. 

Another critical performance criterion of chemical detection systems is detection speed, 
or discovery time. For such a detection system to be practical, it must be able to detect and 
recognize a predetermined quantity of contraband in a short period of time. For example, 
detection of 1 Kg of explosive within about 1 second would practically allow for the use of the 
system in most any application. Electronically, this level of performance requires a high data 
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accumulation rate: at least 10,000 gamma events must be fully measured and processed per 
second (after rejecting associated "noise") in order to have a large enough statistical sample to 
recognize the explosive by gamma ray analysis. Prior art contraband detection systems using 
gamma rays have exhibited detection times on the order of 1 hour for 1 Kg of explosive, or 3,600 
times longer than the desired 1 second previously described, thereby rendering them impractical 
for many applications. This poor detection time is generally caused by two independent factors; 
(i) the "dead time" of the gamma detector; and (ii) the "piling up" of coincidences when more 
than one pulse arrives within the resolution time, which causes accidental coincidences which are 
indistinguishable from the true coincidences. These factors are discussed in greater detail below. 



The "dead time" of an HPGD is the charge collection time, typically 200 ns. This value 
implies an ideal (non-random) maximum counting rate for the detector of 5 million events per 
second (i.e., 1 second/200 E-09 seconds per event = 5 E06 events). However, a practical 
maximum counting rate, or one taking into account randomness, is about 100 times lower; 
roughly 50,000 events per second, assuming the same energy resolution (i.e., 0.1%). Prior art 
detection systems have been restricted by these limits, since the dead time has heretofore been 
considered to be an intrinsic property of the detector. 



"Pile-up", also known as the accidental coincidence rate between two detectors, may be 
represented by the following relationship: 



Dead Time 



Pile-Up 



N a - t r x CRi x CR 2 



(Eqn. 1) 



where: 



= accidental coincidence rate 



t r = event resolving time 

CRi = instantaneous counting rate in detector 1 

CR2 = instantaneous counting rate in detector 2 
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Although Sodium Iodide has good resolving time (on the order of 3 ns), pile-up will begin to 
take place when the product of two counting rates reaches roughly 10 8 ; that is, the random 
coincidence rate will be comparable to the true coincidence rate at that point, and quadratically 
exceed the true rate above that level. It can be readily shown that this corresponds to a neutron 
production rate of 1 E06 neutrons into 4n steradians (i.e., an entire sphere) per second which 
yields about 3,600 times less gamma rays than is required for the desired 1 second contraband 
recognition previously discussed. Hence, a detector system that operates at about 3,600 times 
the rate of prior art systems (while maintaining the accidental count level below a predetermined 
value, such as 10%) would be needed to produce a 1 second recognition time. 

It is further noted that prior art detection systems utilize detector-to-detector coincidence 
circuitry (i.e., an event occurring at one detector is compared to an event occurring at another 
detector), which further exacerbates the aforementioned problems. 

Hence, in sum, no systems or techniques presently in existence provide an effective, 
accurate, and safe method for the non-intrusive detection, location, or analysis of deleterious 
chemical compounds regardless of their physical location or container. A system and method is 
needed by which objects and areas may be rapidly and accurately examined using non-intrusive 
means to determine 1) the three-dimensional (i.e., "X, Y, Z") or spatial location of suspect objects 
contained within, and 2) the chemical identity of the suspect object as an explosive, nerve agent, 
chemical weapon, or item of contraband. 

Summary of the Invention 

The preferred embodiments of this invention relate to an apparatus and method for 
detecting, locating, and analyzing minute or microscopic quantities of chemical elements and 
compounds located within a test subject using subatomic particle activation. 

One aspect of the present invention relates to a system which provides non-invasive 
stoichiometric detection and imaging of chemical elements and compounds in a material to be 
analyzed. The system comprises a particle generator which generates a plurality of first 
subatomic particles and a plurality of second subatomic particles at a target position which is a 
first distance from the material to be analyzed. The system further comprises at least one photon 
detector capable of detecting photons resulting from irradiation of the material by the first 
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subatomic particles and generating a plurality of first electrical signals. The system further 
comprises a particle detector array comprising a plurality of particle detectors, the detector array 
at a second distance from the target position, and the second distance being larger than the first 
distance. The particle detectors are each being capable of detecting at least one second 
subatomic particle from the particle generator and generating a plurality of second electrical 
signals. The system further comprises an analyzer operatively connected to the particle detector 
array and the photon detector. The analyzer comprises a processor that filters the plurality of 
first electrical signals so as to produce a plurality of filtered electrical signals. The analyzer 
further comprises a plurality of electronic coincidence circuits which detect coincidences 
occurring between the plurality of filtered electrical signals and the plurality of second electrical 
signals. 

W Brief Description of the Drawings 

OS Fig. 1 is a side view of a first embodiment of the chemical compound detection and analysis 

S apparatus of the present invention. 

ill Fig. 2a is a perspective view of the apparatus of the present invention adapted for use in 

g examining airport baggage handling cars. 

J;;; Fig. 2b is a perspective view of the apparatus of the present invention adapted for use on a 

M land mobile vehicle. 

S ; Fig. 3 is a perspective view of a first embodiment of an alpha particle detector array 

according to the present invention. 

Fig. 4 is a side view of a second embodiment of the chemical compound detection and 
analysis apparatus of the present invention, utilizing multiple targets and an applied magnetic field. 

Fig. 5 is a side view of a third embodiment of the chemical compound detection and analysis 
apparatus of the present invention, utilizing multiple targets and a single, spatially broadened beam. 

Fig. 6 is a perspective view of a dual-source neutron activation apparatus according to the 
present invention. 

Fig. 7 is a perspective view of a fourth embodiment of the chemical compound detection 
and analysis apparatus of the present invention, as might be attached to a land-mobile vehicle useful 
for land mine or buried explosives detection. 
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Fig. 8 is a block diagram illustrating one embodiment of the gamma spectrum analog-to- 
digital conversion and filtration process of the present invention. 

Fig. 9 is a prompt gamma ray emission spectrum obtained from a typical prior art fast 
neutron activation detection system. 

Fig. 10 is a prompt gamma ray emission spectrum (background subtracted) obtained by 
irradiating a urea test specimen with fast neutrons and detecting prompt gamma emissions using a 
Germanium crystal detector gated at 40 ns. 

Fig. 1 1 is a prompt gamma ray emission spectrum (no background subtraction) obtained by 
irradiating a urea test specimen with fast neutrons and detecting prompt gamma emissions using a 
Germanium crystal detector gated at approximately 6 ns. 

Figs. 12a and 12b are sample prompt gamma ray spectra showing the peak and three 
discrete bandwidths analyzed as part of one embodiment of the chemical identification method of 
the present invention. 

Fig. 13 is a graph illustrating the method of gamma event rise time measurement according 
to the present invention. 

Fig. 14 is a bar graph of a digitized spectrum showing the relative constituent elements and 
concentrations of several common chemical compounds. 

Figs. 15a and 15b are diagrams of the C:H:0 and C:N:0 Dalitz triangles, respectively, 
showing a variety of common chemical compounds. 

N Fig. 16 is a prompt gamma ray spectrum illustrating the chlorine "peaks" associated with the 
interaction of thermalized fast neutrons and chlorine atoms commonly present in cocaine muriate (in 
rice). 

Fig. 17 is a block diagram of the method disclosed in the present invention for evaluating 
and calibrating FNA devices using the discovery time constant Td. 

Fig. 18 is side view of the stoichiometric neutron microscope apparatus with a 
magnification factor i = 380X. 



Detailed Description of the Preferred Embodiments 
Reference is now made to the drawings wherein like numerals refer to like parts throughout. 



Detection, Location, and Analysis Apparatus 

Fig. 1 shows a first embodiment of the particle detection and analysis apparatus of the 
present invention. As shown in the Figure, an accelerated beam 1 0 containing one or more 
subatomic species (here, various ionized isotopes of hydrogen including deuterium and tritium) are 
used to bombard one or more specially constructed targets 18, thereby generating streams of 
subatomic particles (fast neutrons 14 and alpha particles 16) which simultaneously, in pairs, 
emanate from the target(s) 18 in substantially opposite directions. A conventional or advanced 
charged particle accelerator of the typ^ g well known in the art, such jis the Model A-71 1 accelerator 
manufactured by the MF Physics Corporation, is used as the source 15 although a variety of 
different such sources may be used with equal success. In the present embodiment, the source 15 is 
operated in a continuous direct current (i.e., DC) mode such that excitation particles are incident on 
the target(s) continually which may or may not be modulated into "long" (0.1 to 10 sec.) discrete 
time intervals, although other schemes may be used. In the present context, the term "long" is used 
with respect to the coincidence resolving times described later herein, which are on the order of 1 to 
100 nanoseconds. It should also be noted that while neutrons and alpha particles are generated from 
the target(s) in the present embodiment, other subatomic particles or emissions with desirable 
properties may also feasibly be used to produce the desired result within the test subject. The 
object or area being examined 20 is exposed to the generated fast neutron flux 14, the energetic 
neutrons of which interact with the carbon, nitrogen, or oxygen bearing (C:N:0) molecules of any 
explosives, chemical weapons/nerve agents, or contraband 22 within the subject 20 thereby 
generating prompt quanta in the form of gamma rays 24 which are characteristic of C, N, and O and 
are know as "signatures" of the aforementioned substances. It is further noted that the energies of 
these gamma rays 24 have substantially discrete values. One or more high -resolution gamma 
detector s 26 of the t y pe well known in the art are placed relative to the subject being examined 20 to 
detect these emitted quanta, which are subsequently analyzed to identify the molecules resident in 
the test subject. Additionally, one or more scintill ation (or comparable )_deLectors 28 of the type 
well kmvj^m^e^art are placed within the solid angle of the alpha particle flux 16 in order to detect 
alpha particles in coincidence with the gamma events detected by the aforementioned gamma 
detectors 26, thereby permitting two separate determinations: 1) that of the gamma-alpha 
simultaneous timing (rejecting non-synchronous events), and 2) spatial positioning. High purity 
(80%) germanium detectors of the N-type (neutron resistant) produced by ORTEC Corporation 
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are used as the gamma detectors 26 in the present embodiment. Note that while Germanium 
crystal detectors are preferred (as further described below), it can be appreciated that other types of 
high-resolution detectors, for example those utilizing Xenon, may be used to accomplish the desired 
functionality of coincident alpha and prompt gamma detection for both the timing and spatial 
location of the chemical compound within the test subject. 

The targets 18 of the present embodiment are constructed of a scandium titride layer 
deposited on a copper (Cu) substrate, both substances which are well known in the material 
sciences, the targets having the desirable property of generating a stream of neutrons (neutrally 
charged nucleons) and alpha particles (ionized helium nuclei) when properly excited by the incident 
deuteron/tritium ion beam 10. It can be appreciated, however, that other types of targets and 
materials may be used in this application. 

Neutrons are created in the target(s) according to the following exemplary reactions: 

!H 2 + iH 3 6 2 He 4 + 0 n ! + 14 MeV 

iH 2 + jH 2 6 2 He 3 + on 1 + 3.27 MeV 

]H 3 + iH 3 6 2 He 4 + 0 n J + 0 n ] + 1 1 33 MeV 
Note that the energy of the incident deuteron/tritium ion beam 10 must be sufficient to overcome 
any coulombic interaction with the positively charged nuclei of the target material atoms. 
Deuteron/tritium ion energies of 0.05 MeV or greater have been found sufficient for this purpose. 

Fast neutrons having energies on the order of 14 Mev are utilized to bombard the test subject 
18 in the present embodiment due to their desirable scattering properties (i.e., inelastic scattering 
with nuclei) and ability to penetrate significant thicknesses of common substances such as steel, 
soil, sand, lead, earth, and slabs of water up to approximately 50 cm in thickness (1/e interaction 
length). Note that the cross section (in mbarns) for gamma production in C, N, and O by 14 MeV 
neutrons is nearly independent of the neutron energy at that energy level; thus, the relative 
concentrations of these elements can be obtained to a high degree of accuracy without knowing the 
actual collision energy. This is in contrast to the lower neutron energies of many prior art systems, 
which have cross sections which vary much more significantly with neutron energy, thereby making 
it practically impossible to calculate the relative chemical contributions without knowledge of the 
precise collision energy. Despite these considerations, however, it will be appreciated that neutron 
energy levels other than 14 MeV (and even multiple energy levels) may be used based in the present 
invention on the desired system operating characteristics. 



The neutron and alpha particle beams 14, 16 released by the target are distributed spatially 
throughout a given solid angle <|> (measured in steradians) which is related to the angle of incidence 
of the charged particle beam to the target, the axis of each beam (corresponding to the highest 
neutron or alpha particle flux) being substantially co-linear with the other, yet opposite in direction 
relative to the target 18 (see Fig. 1). Each target may be either fixed or independently steerable 
(adjustable) in relation to the excitation beam 10 and test subject 20 via a conventional electro- 
mechanical positioning device 21, although it can be appreciated that any variety of arrangements 
may be used. Such positioning devices may be manually controlled, or alternatively automatically 
controlled via inputs such as those from the signal processing components described later herein. In 
this fashion, the resulting neutron and alpha particle beams 14, 30 may be adjusted to provide the 
desired neutron/alpha flux in a given solid angle, such as for scanning purposes. 

In the present embodiment, alpha particle detection is accomplished via an array 34 of 



scintillation or comparable detector elements 36 having a'known spatial relationship to each other. 



Fig. 3 shows an exemplary configuration consisting of a square array (8x8) of 64 discrete detector 
elements. Each of the 64 detector elements is viewed by one photomultiplier ("PM") tube or by one 
of the 64 "segmented anodes" built into one large PM tube, each anode acting as an independent 
counter. This array 34 is placed in a known position relative to the target(s) such that each detector 
element subtends a given solid angle § with respect to its associated target 18. Accordingly, based 
on the substantially co-linear relationship between the fast neutrons and alpha particles emitted from 
the target upon excitation, individual alpha detection events can be directly correlated with the locii 
of the neutron-induced events occurring in the test subject within a similar but inverse solid angle 
(see Fig. 1). 



It should also be noted that in the present embodiment, borated polyethylene elements 38 (or 
those constructed of comparable neutron absorbing or moderating material) are used to shield 
personnel and equipment adjacent to the apparatus from the deuteron/neutron/alpha radiation 
generated within the system, and further to collimate the neutron beam 14 generated by the target(s) 
if desired as shown in Fig. 1. 

Note that the present embodiment (as well as those later discussed) may be adapted to a 
variety of different applications and geometries including, inter alia, land mine detection and 
identification, artillery shell analysis, or as shown in Fig. 2a, analysis of airport baggage carriers for 
contraband. In the context of detecting buried or otherwise shielded explosives such as land mines 
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or artillery shells, the present invention may be readily adapted to a land mobile vehicle of the type 
well known in the mechanical arts (Fig. 2b) which is either controlled directly by an operator, or 
remotely controlled in order to facilitate surveying large portions of terrain with relative efficiency. 
It will be recognized that the greatly reduced detection time associated with the present invention 
(on the order of a few seconds) makes such land mobile applications feasible, whereas such 
applications were not feasible using prior art particle activation systems due to their comparatively 
long detection times (characteristically on the order of tens of minutes or more). Such a land 
mobile vehicle would ideally have the source(s) 15, target(s) 18, and detectors 26, 36 physically 
removed from the chassis of the vehicle, such as on an extensive boom or arm (not shown), so that 
sufficient distance between these components and the vehicle could be maintained and the chance of 
damage to the vehicle reduced, although other configurations are possible and within the scope of 
the present invention. 

In a second embodiment of the aforementioned apparatus (shown in Fig. 4), individual 
atomic species resident in the excitation beam 10 are separated through the application of a 
magnetic field 40 induced along the beam path in order to permit the excitation of more than one 
target 18. Specifically, the excitation beam produced by the aforementioned source 15 contains a 
plurality of atomic species including deuterons and tritium ions, each having different atomic mass 
number. As is well understood in the physical sciences, a charged particle passing through a 
magnetic field experiences a deflecting force, the magnitude and direction of which is determined 
by the particle's charge and mass, and the strength and direction of the magnetic field vector at that 
given location. Hence, particles of different atomic mass but of the same kinetic energy can be 
deflected along curved paths of different radii utilizing the same magnetic field. In the present 
embodiment, two atomic species are deflected using a perpendicular magnetic field 40 in such a 
manner as to permit impact of the different species on two different targets 18. The magnetic field 
is generated through use of a conventional or superconducting electromagnet 42 of the type well 
know in the art which is placed in direct proximity to the excitation beam path 10, although other 
methods and arrangements (such as electrostatic depletion) may be utilized with equal success. The 
charged excitation particles strike two spatially disparate targets of the same construction as 
previously described, thereby permitting both resulting fast neutron beams 14 to be directed toward 
the same test object or area simultaneously. This "binocular" apparatus permits enhanced spatial 
resolution of the gamma-emitting chemical compound within the test subject via triangulation; i.e., 
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the emission location is fixed by computing the intersection of two lines (or subtended solid angles) 
within the test subject 20. This way, the usual determination of the Z-axis coordinate by measuring 
time-of-flight of each neutron is avoided and precise X,Y,Z values are obtained directly. Spatial 
resolution on the order of 1 cm or less is practically achievable using, inter alia, the binocular 
apparatus of the present invention. 

It should be noted that alternatively, the excitation beam 10 may be split into multiple beams 
and directed into multiple targets using (i) an electrostatic septum, as is well known in the field of 
particle physics, or (ii) by broadening the beam spatially and utilizing multiple targets to intersect a 
fraction of the beam, as shown in Fig. 5. Using these approaches, the beam 10 may be comprised of 
one or multiple species. 

In yet another embodiment of the invention shown in Fig. 6, two or more deuteron sources 
15 are utilized to excite a plurality of targets 18, thereby inducing the emission of neutrons for the 
targets. This embodiment obviates the need for the application of the aforementioned magnetic 
field 40 (and the supporting structure necessary to generate such field) or broad particle beam. The 
multiple deuteron sources 1 5 are physically placed so as to optimize the ability of the device to 
spatially locate the chemical compound(s) of interest within the test subject 20, and are similarly 
operated in either a continuous ("long") or modulated mode (such modulated pulses being 
substantially coincident or staggered in time, depending on the desired properties of the system). 

Referring now to Fig. 7, another embodiment of the chemical compound detection, 
location, and analysis apparatus of the present invention is shown. A plurality of gamma 
detectors 26 are placed in known physical proximity to the test subject, and one or more alpha 
particle detectors (not shown) are placed in general proximity to the target(s). The target is 
excited by a DC deuteron pulse of known duration (typically 0.1 to 20 sec), and the timing 
difference between a given alpha particle detection event (or series of events) as detected by the 
alpha detectors 28 and the induced prompt gamma event(s) detected by the gamma detectors 26 
is measured to provide a coarse determination of the axial (i.e., line-of sight, or "Z-axis") 
distance 44 between the activated chemical compound and the target 18. Such determination is 
performed by a digital signal processor (or microprocessor) running an algorithm in the present 
embodiment; although other arrangements may be used. Known neutron and alpha particle 
velocities, electronically-induced detection delays, prompt gamma emission delays, and 
excitation beam pulse duration are considered in order to produce this coarse estimate. Further 
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refinement of the coarse axial position previously calculated, as well as an off-axis (i.e., "X-Y") 
position determination, are provided through subsequent multi-node geometric analysis (using 
any number of lines of position, such as triangulation) of data provided by the aforementioned 
gamma detector array 26. Specifically, by knowing the angular and spatial relationships of the 
gamma detectors to each other and the target, the location of a prompt gamma emitter (i.e., 
chemical compound) can be determined through analysis of the counting rate and/or timing 
differences between individual detection events. This arrangement also has the advantages of 
not having to measure neutron time-of-flight (TOF) or neutron emission/incidence angle. 



Coincidence Detection 

Referring now to Figs. 3 and 8, an improved coincidence detection apparatus and method is 
described. An array of alpha particle detectors such as that shown in Fig. 3 and previously 
described herein is formed in general proximity to the neutron-generating target(s) previously 
described. Pre-selected individual gamma spectral lines (corresponding to the emission of prompt 
gammas from a given type of atom within the test subject) are electronically correlated with each 
alpha detector array element 36 by the analyzer 39 to form a plurality of parallel coincidence 
circuits. Specifically, gamma detection events remaining in the post-filtration gamma spectrum 
(Fig. 8) are correlated with alpha detection events for each alpha detector array element using a 



nominal 10 msec, gating window. In on ^ embodiment , approximately twenty (20) different discrete 
gamma lines from one gamma counter 37 and detector 26 (corresponding to the emission lines of 
the elements of concern) are selected and electronically placed in coincidence with signals from 
each of the 64 elements of the alpha detector array 34. This arrangement produces a large number 
of parallel coincidence channels. For example, in the present embodiment, a set of four gamma 
detectors 26 and 64 alpha detector elements 36 are used. As will be described in further detail 
below, each of the 20 selected spectral lines is decomposed into three roughly 4 KeV bandwidth 
components, thereby creating 20x3=60 independent channels per gamma detector/counter. Since 
four detectors are used, 60x4=240 independent channels are associated with the gamma detector 
array as a whole. When correlated with the 64 alpha elements 36 in the alpha detector array 34, a 
total of 15,360 independent coincidence channels result. Contrast this with a prior art "detector- 
detector" coincidence system, which using the same 4/64 arrays, would produce only 256 
independent channels. This highly parallel coincidence circuit structure of the present invention 
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permits the processing of a substantially increased number of prompt gamma events, since only 
those events correlating to the desired twenty or so spectral lines (C:N:0 in the present 
embodiment) need be correlated to events detected by the alpha array 34, and such events are 
processed with a high degree of parallelism. Accordingly, a higher incident neutron flux 14 (on the 
order of 10 9 to 10 11 n/s-4pi) may be used as compared to prior art systems using detector-detector 
coincidence (and maximum neutron fluxes on the order of 10 6 n/s-4pi), and significantly less time is 
ultimately required to identity and analyze a given compound under the same test conditions. 

It will be recognized that precise gamma energy determination is important in the present 
invention in order to identify the parent atom of gamma rays chosen for analysis. Several factors 
influence the selection of gamma ray peaks and spectral lines used for the analysis. These 
factors include (i) gamma energy; (ii)cross section; (iii) cascade vs. photo-peak; (iv) proximity 
and overlap; and (v) single/double escape peaks. These factors are discussed in greater detail in 
Appendix A attached hereto. It will also be appreciated that the coincidence and detection circuitry 
of the present invention may optionally utilize corrections for noise and detector efficiency, one 
approach for described in ^Appendices B and C, respectively. 

In a sec^nd^en^odimenTo^the coincidence detection apparatus of the present invention, 

r ■ i,, „||. i,. .il - l l-l 

38,400 coin^*dence-circutts"(all accommodated by the signal processing capability of one signal 
processing IC of the type well known in the electronic arts) are formed. Specifically, the alpha 
detector array 34 is segregated into 640 discrete detector elements 36, and the output of the 
gamma detector(s) 26 is separated into 30 discrete energy bandwidths. Separation of the gamma 
detector output into bandwidths as described results in each gamma detector bandwidth acting as 
a separate detector electronically. Accordingly, with the 640-element alpha array 34 and a single 
physical gamma detector (30 "electronic" gamma detectors), 640x30 = 19,200 coincidence 
circuits may be formed. Lookup tables are used in conjunction with a standard random access 
memory (RAM) within the analyzer 39 of the present embodiment to facilitate rapid processing 
of this large number of coincidence circuits. In addition, each coincidence circuit in the present 
embodiment is provided a parallel coincidence circuit with an artificially imposed delay to 
account for the aforementioned "accidental" coincidences. This amounts to an additional 19,200 
coincidence circuits. Hence, the total number of coincidence circuits used in this second 
embodiment is 38,400 (19,200 + 19,200). 
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Chemical Identification 

Referring now to Figs. 9 through 14, an improved gamma ray detection, filtration, and 
analysis apparatus and method is described. As shown in the Figures, high purity Germanium (Ge) 
crystal detectors (HPGDs) 26 are used to detect prompt gamma ray emitted from within the test 
subject. The Germanium detectors provide enhanced gamma energy resolution (on the order of 
0.1%) unlike the more common Sodium Iodide detectors (6-10%) typically used in prior art 
systems, thereby allowing discrimination of the multiple C:N:0 or other spectral "lines". Fig. 9 
illustrates a typical gamma spectrum obtained using a prior art Nal detector. Note the breadth of the 
detected peaks, which is due to the comparatively poor energy resolution of this type of detector. In 
contrast, the high-resolution gamma spectrum 50 produced from the Ge detectors of the present 
invention provide greatly enhanced energy resolution, as evidenced by the narrow spectral peaks or 
lines illustrated in Figs. 10 and 11. In the present embodiment^tl^gamm^pectrum 50 is 
electronically converted to a digital representation (Fig. 14) using a conventionai^analog-to-digital 
(A/D) converter 3 1 of the type well know in the electronic arts. Each spectraHi jie 52 is as signed a 
discrete binary value ("bin") representing its gamma energy level.<^Kmwn_spectral lines (bins) 
associated with carbon, nitrogen, and oxygen for the chosen type/energy of incident particle stream 
are then identified as described in Appendix A and further prpcessed~while_other unrelated lines 54 
are electronically filtered (using, for example, ^co^^ To digitize the 

gamma spectral lines, three (3) coincidence circuits are used per each line. Specifically, in one 
\ embodiment of the present invention, amplitudes of each of the aforementioned approximately 20 
spectral lines associated with C:N:0 are determined through analysis of three equal 4 KeV 
bandwidths within the spectral line; namely, one "peak" bandwidth and two "shoulder" bandwidths 
(see Figs. 12a and 12b). The peak-to-shoulder difference(s) are used to determine the amplitude of 
the peak for purposes of further analysis in the Dalitz plot. As previously discussed, Germanium 
detectors characteristically have a slower response time than other types of detectors (such as 
Sodium Iodide crystal), thereby having a correspondingly lower temporal resolution. It is assumed 
that HPGDs can process a maximum event rate (including random events) on the order of 50,000 
counts/sec. This slower response rate is compensated for in the present invention through the use of 
electronic processing of signals from the detector which effectively varies the response time 
constant to a lower value. Specifically, in the present embodiment, that portion of the HPGD 
signal corresponding to a fraction of the rise time of the gamma event is used to determine the 
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time resolution. This rise time is typically in the range of 1.5 to 4 ns, and is measured from a 
point 10% above the baseline prior to the event to a point 10% below the peak value of the event, 
as shown in Fig. 13. The rise time signal processing is accomplished . via a constant fraction 
(7 discriminator (CFD) which is ^e^kno^n-in the signal processing and nuclear detection arts. 



Charge collection in the present embodiment is further stopped electronically ("gated") at 20 ns. 
Using this arrangement, the effective maximum count rate of the Ge detector is substantially 
increased, since the detector "dead time" is reduced, and temporal resolution increased. 

Figs. 10 and 1 1 are exemplary gamma spectra which depict the effect of gating on spectra 
obtained from a given sample of a chemical compound. Fig. 10 illustrates a gamma spectrum 
obtained from a urea specimen using a nominal 40 ns gating interval, and background 
subtraction. Background subtraction is performed in the present embodiment using the method 
described in Appendix D hereto, with reference to Fig. 12b. 

Fig. 11 illustrates a gamma spectrum of the same urea specimen taken using a 5.9 ns 
gating interval, with no background subtraction (background subtraction is precluded in the 
spectrum of Fig. 11 due to the narrow gate). When comparing the spectra of Figs. 10 and 11, it 
is readily noted that the signal-to-noise ratio (SNR) is significantly enhanced in the spectra of 
Fig. 11, due primarily to the increased temporal resolution afforded by the shorter gating 
interval. Hence, an appreciable increase in performance is obtained by using a shorter gating 
interval, without complex noise suppression or software manipulation. 

Experimental data obtained by the applicant herein indicates that counting rate increases on 
the order of 100% to 400% over non-gated HPGDs are possible when using the techniques 
described herein. Furthermore, the HPGD of the present invention can meet both requirements 
for practical contraband detection and identification simultaneously; i.e., high energy resolution 
and high temporal resolution. Note that while the use of electronic gating at 20 ns herein results 
in a decrease in energy resolution (i.e., from about 0.1% without gating to about 0.3% with 
gating), the resultant energy resolution is more than sufficient for the purposes of contraband 
identification according to the method of the present invention. It will further be recognized that 
while a CFD performing electronic gating at a nominal interval of 20 ns is used in the present 
embodiment, other gating intervals, types of circuits, and techniques may be employed to 
measure and utilize the desired portion of the gamma detector signal. 
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It is further noted that by using the foregoing filtering and gating techniques, the 
computational burden on attached signal or data processing equipment is greatly reduced, and more 
individual gamma detection events can be processed per unit time, since only relevant C:N:0 (or 
other) spectral lines survive the filtration stage and need be, correlated with detected alpha particle 
events/ Processing of the resulting spectral signals 56 in real time may be accomplished using any 



variety of conventional digital slgfral processing devices such as a "Lookup Table" LeCroy Model 
2373 or comparablFwhTch~pTovides a high data processing rate. 

Appendix E hereto describes the aforementioned gamma detection, selection, and 



Dalitz Triangle 

Referring now to Figs. 15a and 15b, an improved method for identifying chemical 
compounds using particle activation is described. The filtered digital gamma detector output signal 
described in the preceding paragraph is input to a signal processing apparatus (such as the 
aforementioned digital signal processor) running. an improved identification/classification algorithm 
that reduces a three-dimensional problem to a two-dimensional representation. This algorithm is 
based on the well known "Dalitz Triangle" 60, 62, which correlates the concentrations of three 
elements within a compound. See Appendix F. Specifically, the length of the normal to each side 
of the equilateral triangle 60, 62 is proportional to the square of the atomic density of each of the 
three elements. The algorithm of the present invention takes the filtered digital gamma spectrum 56 
which is collected over a given integration interval, computes the relative proportions of the 
constituent elements based on the binary representations of amplitude of the spectral lines (which 
correlates to the number of prompt gamma events detected), and mathematically generates a vector 
64 normal to the appropriate side of the triangle for each element being analyzed. See Fig. 15b. 
The point of intersection 66 (or closest point of approach) for all three normal vectors is calculated 
using well known geometric techniques, thereby defining a single point or range. This point/range 
is then compared to a "library" of points/ranges associated with various known chemical 
compounds of interest (such as may be stored in tabular form within a conventional random access 
memory (RAM)) to identify the compound detected. Subsequent spectra collected from the gamma 
detectors over a given sample period are then used to statistically increase the confidence level of 
the identification result. Statistical error bands 68 (such as a given number of standard deviations on 




processing utilized by the present invention in additional detail. 
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a normal or Gaussian distribution) may be used to establish confidence criteria. Such data may also 
be displayed visually to the system operator, thereby providing him a visual indication of the 
confidence of a given analysis. For example, the Dalitz triangle may be displayed in graphical 
form, with the library of known chemical compounds displayed as discrete points within the triangle 
60,62. As successive spectra are obtained and analyzed by the system during its integration time, 
they may displayed as other discrete points on the triangle based on their relative proximity to the 
library point (compound) of concern. It can be appreciated that a large variety of different display 
formats (including tables, bar graphs, scatter plots, and the like) may be utilized to convey the 
desired information, the generation and rendering of all of such formats being well known within 
the relevant arts. 

Auto-Thermal Neutron Activation 

An apparatus and method for the simultaneous or sequential application of both thermal and 
fast neutrons to identify a chemical compound is now described. For the purposes of the present 
discussion, thermal neutrons are considered to be those neutrons with a total kinetic energy level 
substantially less than those of fast neutrons. For example, thermal neutrons used in the present 
invention may have energies on the order of 0.025 eV, while fast neutrons may have energies on the 
order of 14MeV as previously described. It will be appreciated that neutrons of multiple energy 
levels may be present within the object under examination when irradiated by fast neutrons, and that 
the detection and analysis of quanta emitted as a result of these various neutrons is within the scope 
of the present invention. 

Referring now to Fig. 16, an exemplary gamma ray energy spectrum associated with the 
cocaine molecule (C17H21NO4.HCI, also known as cocaine muriate) is shown. Cocaine muriate is a 
form of cocaine commonly used in the illicit drug trade. The spectrum of Fig. 16 is produced 
without the coincidence requirement by detecting both prompt gamma rays resulting from the 
irradiation of the cocaine muriate molecule with fast neutrons as previously described, and the 
delayed gamma rays from the capture of thermal neutrons in the nuclei. Generally, two primary 
mechanisms are involved in producing a spectrum such as that of Fig. 16: (i) fast neutron activation 
of, and subsequent prompt gamma emission from, the carbon, nitrogen, and oxygen atoms within 
the illicit substance (cocaine muriate in this example); and (ii) thermal neutron activation of other 
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"pointer" atoms present in the illicit substance (here, chlorine), also known as "auto- 
thermalization' ' . 

Using the apparatus of the present invention described with reference to Fig 1 above (or 
other embodiments as described herein), the object suspected of containing illicit substances is 
irradiated with fast neutrons generated by a neutron source. These incident fast neutrons are 
moderated or thermalized within the material of the container, within any material surrounding the 
illicit substance, and by the illicit substance itself. For example, the aforementioned cocaine is often 
times packed in ordinary rice as a means of confounding detection systems; when an object 
containing cocaine and sugar is irradiated with fast neutrons, the hydrogen and carbon atoms present 
in the cocaine and rice slow the fast neutrons to thermal energy levels via multiple inelastic 
scattering events. The thermalized neutrons subsequently interact with the chlorine atoms in the 
cocaine, these "pointer" atoms having a comparatively large thermal neutron scattering cross- 
section of about 33 barns, and produce a series of tell-tale delayed gamma emissions having energy 
levels centered at approximately 2.69, 6.62, 6.86, and 7.42 MeV as shown in Fig. 16. Note also that 
in addition to rice, common drug packing substances such as sugar (sucrose) or coffee which are 
also rich in hydrogen and carbon atoms are also good "auto-moderators" of fast neutrons. 

In the apparatus and method of the present invention, the aforementioned tell-tale gamma 
emissions are detected and counted using the HPGDs and associated scaling circuitry described 
above, and when a desired statistical confidence level met (e.g., 500 counting events in the selected 
peak), a signal or "flag" generated to indicate the possible presence of cocaine muriate. The 
chlorine peak with the highest signal to background ratio (e.g., that occurring at 7.4 MeV) is 
selected as the flag in the present embodiment, although it will be appreciated that other peaks may 
be used. Note that the signal to background ratios for the 6.86 and 7.42 MeV peaks are roughly an 
order of magnitude higher that the ratios obtained from conventional prior art thermal neutron 
activation systems under the same circumstances (e.g., 2:1 to 3:1, versus 0.2:1 to 0.3:1 for the prior 
art systems). This roughly ten-fold increase in signal-to-background in the present invention is 
attributable to the auto-thermalization of fast neutrons within the contraband and its packing 
material. 

In parallel with the foregoing analysis of the gamma spectrum attributable to the chlorine 
atom, the prompt gamma spectrum resulting from the fast neutron scattering is analyzed as 
previously described herein. Hence, the gamma spectrum resulting from fast neutron activation is 
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used to refine or confirm the warning signal provided by the thermal neutron induced gamma 
spectrum. Such confirmation is needed, since arguably many materials of a non-illicit nature 
contain chlorine atoms as well. Note that no gamma/alpha coincidence or spatial analysis of the 
thermal neutron induced gamma spectrum is performed in the present embodiment; this facilitates 
very rapid and simple processing, and allows the warning signal to be generated before the 
completion of the empirical analysis of the fast neutron induced quanta. 

It will be recognized that while the foregoing apparatus and method is described with 
respect to cocaine muriate and its chlorine atoms, other illicit substances may be analyzed and 
detected using the present invention. For example, heroin (diacetylmorphine hydrochloride 
monohydrate, chemical formula C21H23NO5.HCl.H2O) includes chlorine atoms which may be used 
to generate a gamma peak indicative of the possible presence of heroin. Other "pointer" atoms 
having a substantial thermal neutron cross-section may conceivably be used as well, such as 
potassium, titanium, vanadium, chromium, manganese, cobalt, silver, cadmium, indium, thallium, 
tungsten, or mercury. Additionally, while the analysis of the gamma rays resulting from fast 
neutron activation is preferably performed in parallel or contemporaneously with the thermal 
neutron gamma analysis, it will be appreciated that such analysis my be sequential or serial in 
nature, or even performed before the thermal neutron analysis if desired. 

It will further be recognized that the foregoing apparatus and method has the advantage of 
analyzing quanta resulting from thermal neutrons generated within the object under examination, as 
opposed to predominantly incident thermal neutrons as used in prior art neutron activation systems. 
In these prior art systems, fast neutrons are typically made to hit a moderator (usually water or 
paraffin) near the target where they are thermalized and lose directional orientation. As a result, the 
thermal neutron flux emitted from the moderator toward the object under examination is only a 
small fraction of the fast neutron flux incident on the moderator. This effect greatly reduces the net 
gamma flux out of the object being examined, thereby reducing counting rate and system efficiency. 
Note that in such prior art systems, a small fraction of fast neutrons incident on the moderator may 
escape the moderator without being thermalized, and therefore may go on to be thermalized within 
the object under examination. This is in contrast to the invention disclosed herein, in which 
effectively all thermalizations occur within the object under examination. 

Additionally, the penetration depth of the fast neutrons used in the present invention is 
greatly enhanced (on the order of 1 meter), as compared to the typical penetration depth (roughly 2- 
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3 cm) for incident thermal neutrons. This increased penetration depth is important to the efficacy of 
the system, since most contraband is secreted within other objects which often times have 
significant mass and density, thereby significantly attenuating thermal neutrons within a very short 
distance long before they hit the contraband. Furthermore, because only a small fraction of incident 
thermal neutrons might reach the hidden contraband when using one of the aforementioned prior art 
systems, the resulting gamma flux out of the object under examination is accordingly small, thereby 
necessitating either a very high incident thermal neutron flux, or a very long counting/integration 
time. 

System Calibration and "Figure of Merit" 

An improved method for measuring the efficacy ("figure of merit") of a particle activation- 
based detection/analysis system and calibrating same is now described. Referring now to Fig. 17, a 
discovery time constant (T d ) is mathematically defined as follows: 

M 

Where: 

Td = Disc. Time Constant (s) 

Li = Target - Chemical Dist. (m) 

L2 = 2 Detect - Chemical Dist (m) 

M = Mass of Chemical (Kg) 
Note that the quantities LI and L2 are squared due to the solid angles subtended by the system 
detectors. This time constant is a measure of the time required to perform a chemically specific 
identification of a certain mass of chemical compound at a certain distance from the particle source 
(target), and certain distance between the chemical compound and gamma detector(s), with a 
prescribed statistical confidence level about the relevant point on the aforementioned Dalitz triangle 
60,62. In practice of this method 700, the relevant distances to the chemical compound and 
detector(s) are measured or calculated, and the mass M determined as shown in step 702 of Fig. 17. 
Next, the desired confidence level is selected in step 704. This confidence level is mathematically 
determined through statistical analysis of the spectral data produced on successive operation of the 
system during a given testing interval, as previously described. The known mass of a given 
chemical compound is then analyzed in step 706, thereby experimentally determining Td for this set 
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of constraints. Next, in step 708, the time constant is optionally normalized. For example, the 
system may require 5 seconds to identify a 1 kg sample compound to a 95% level of confidence at a 
detector/sample distance of 1 meter in dry air at STP; this nominal or "baseline" value may be 
normalized to 1 second if desired. In steps 710 and 712, the mass and/or environmental test factors 
(such as location of the mass, relative humidity, etc.) are varied and the discovery time again 
measured under these new conditions. These steps 710, 712 may be repeated a number of times 
when, for example, incrementally varying one parameter, or varying many parameters at once, as 
described further below. Next, correction or calibration factors are calculated for the system in step 
714. For example, if the system was normalized using a 1 kg mass in dry air at STP, operation of 
the system in humid air at lower temperature (assuming comparable mass and detector geometry) 
would produce a higher value of T d . This higher value of Td is then ratioed with the nominal Td to 
produce a correction factor for operation of the system in humid/low temperature environments. 
Tables or curves may also be constructed detailing the response of the system as a function of 
varying system parameters and test conditions, or salient combinations of parameters. Lastly, other 
critical system constants, including 1) the rate of system "false positives", and 2) the rate of "false 
negatives" (e.g., missed detections/identifications where an actual compound of interest was 
present) are experimentally determined and correlated to the confidence level previously described 
in step 716. For example, if the false positive rate increases substantially below a confidence level 
of 90%, this value is used as a minimum threshold for system operation. The T d value necessary to 
achieve this confidence level or higher is then specified as the minimum operating time for the 
system when examining objects or areas. 

In addition to the foregoing uses relating to calibration and field operation, experimentally 
generated values of Td for one system (under a certain range of specified test conditions) may also 
be compared to those generated for another system in order to determine the relative efficacy (merit) 
of the two systems at chemical identification and location. This method is very useful in 
standardizing the comparison of two systems of differing operating principle; at present, no system 
of comparison exists, hence manufacturer's claims of performance can not be directly compared. 
Such "figures of merit" may be developed as a function of various critical parameters (such as 
ambient temperature, humidity, interposed shielding, etc.) to assist in determining the best 
instrument/technique (or combination thereof) to apply in a given testing or field application. 
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For brevity, the invention described above is referred to below as an 'atometer' because it 
determines the number of atoms of chemical elements contained in each volume element of the 
substance investigated. The atometer invention renders it possible to non-invasively 
stoichiometrically identify and locate, from a distance, explosive and drug through thick layers 
of steel, soil, rice and other materials, by means of stoichiometric analysis and imaging by 
treating prompt gamma rays produced by 'tagged' fast neutrons. A neutron is 'tagged' by the 
light flash generated by the alpha particle in a scintillation detector, said alpha particle that is 
produced in the same nuclear reaction simultaneously with the neutron and whose direction of 
motion is exactly 180 deg (in center-of-mass) to the neutron direction. Fast neutrons, 
inelastically colliding with nuclei inside the object investigated produce prompt gamma rays 
whose energy is characteristic of the element with which neutron collided. Only the gammas 
that are produced at the "same time" with an alpha are accepted as true events, here, 'same time' 
means within a certain prescribed number of nanoseconds; such timing restriction is 
accomplished by requiring strict coincidence between gamma ray detected in a solid state (high 
energy resolution) gamma detector and the alpha flash in a scintillation detector. An electronic 
system - a parallel processor using multitude of fast coincidence circuits, filters out all gamma 
spectrum except the three lines corresponding to the elements of the explosives: C, N and O; 
digitizes their amplitudes and computes a, b and c in the empirical chemical formula in the form 
C a N b O c ; that is, it deciphers on line the empirical chemical formula of the neutron irradiated 
substance. The imaging is obtained simultaneously by taking advantage of the fact that every 
alpha is emitted at exactly 180 deg (in center-of-mass) to the neutron direction; hence, the x, y 
position of alpha particle in a position sensitive alpha detector determines the neutron direction; 
and the neutron time-of- flight between alpha and gamma provides the z coordinate. Distribution 
of alpha particle counts as a function of x and y in said position sensitive alpha detector provides 
2 dimensional images. 

A critical parameter of atometer is the imaging factor, i, which is the ratio i = L/L 2 , 
whereas L is the distance between the source of alpha particles and detector or array of detectors 
capable of detecting alpha particles, or 'detector arm length'; and L 2 is the distance between 
source of neutrons and object investigated, or 'object arm length', and whereas the alpha particle 
source and neutron source are one and the same source. 

Certain chemical inspection of inaccessible substances require that the image of the 
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substance investigated be magnified by a significant factor in order that its detailed features be 
discerned. For example, in semiconductor industry, the wafer inspection during and after the 
manufacture is required to quantitatively establish impurities, inclusions, and, generally, the 
amounts of micro quantities of the desirable (doping) and undesirable chemicals. Such chemical 
imaging is not feasible with the X ray devices currently used, because X rays are chemically blind; 
they respond only to density changes. Another example is found in biotechnology where the 
analysis of micro quantities of extraneous chemical elements and compounds is essential (in blood, 
lab samples of bacteria and other living organisms or DNA, RNA etc). An important example is the 
need for identification and fingerprinting of the gems such as diamonds, sapphires, emeralds etc. 
which have different inclusions of metals depending on the their provenance. Typical microscopic 
impurities are: Ca, Ti, Cr, Mn, Zn, Sn, Sb, Mg, Al, Fe and Ni. In ruby and sapphire, the impurities 
are greater and are dominated by Va, Ti, Mo, Cr and Fe. Detection and identification of impurities 
in such materials require magnification of the chemical images from i = 1.0X to i = 10,000 X. 

All the current techniques use ionizing particles which can potentially change the crystal 
structure and color. The laser techniques are unreliable when more than one impurity is present and 
slow; increase of the examination speed would require an increase in laser power that can damage 
the gem. Electron microscope cannot penetrate more than about 1 to 10 microns in thickness. In 
contrast, fast neutrons do not ionize the materials in their passage, are not stopped by solids or 
liquids, can penetrate up to 0.5 m and provide stoichiometric information deeply inside the gems, 
and would not change the crystal structure. 

Another need for magnification is the detection of microscopic cracks in metals and other 
materials by fast neutron imaging as neutron radiography. 

The preferred embodiments of the present invention satisfy the aforementioned needs by 
adding an image magnification feature to the atometer invention for detection, location and 
stoichiometric analysis of micro quantities of chemical compounds with the purpose of chemical 
identification of intrusions, doping, inclusions as well as the cracks in metals and composite 
materials undetectable without magnification. 

Referring to Figure 18, magnification is accomplished by making the detector arm Li 
longer than the object arm L 2 so that the imaging factor, i = Li/L 2 becomes larger than unity. 
With i > 1 , the image will be magnified by the same ratio. 
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Figure 18 shows an embodiment for stoichiometric detection inside of a diamond sample 
800 by fast neutron activation. A deuterium beam 810, D + , from ion source 820, is accelerated 
and is incident onto tritiated target 830. Alpha particles 840 are emitted upwards to alpha 
window 850, where they are detected by an array of detectors 860. Fast neutrons are produced at 
the target 830 and emitted downwards towards the object 800 investigated - a diamond (black 
circle) at a distance L 2 . Gamma rays 870 produced by the neutrons in the diamond 800 are 
detected by a High Purity High Resolution Germanium Crystal Detector 880. In Figure 18, the 
object diamond 800 is drawn outside target box 835 for clarity; in practice, it is placed inside the 
target box 835 in which case L 2 = 0.3 cm. Detector distance is Li = 1 15 cm and the image ratio 
is i - 1 15 cm/0.3 mm = 380. Hence, the image of object will be magnified 380 times at the alpha 
window position. Each linear dimension of the object, x 2 or y 2 , at distance L 2 will be magnified 
by the factor 380, thus the image dimension will be x = ix 2 and y = iy 2 . 

The magnification i can be augmented by a straightforward mechanical increase of Li to 
1000X. An additional increase of i is achievable by magnetically lengthening alpha particle 
orbits in a magnetic field produced by a solenoid, the axis of which is at about 80 deg to the 
vertical axes in Figure 18. The field shape can be either non-focusing or focusing. The central 
orbits, known as migma orbits, will be formed as described in Phys. Rev. Lett. 54, 746 (1985), 
which is incorporated by reference in its entirety. Depending on the filed topology, their length 
can be in kilometers. The maximum radius of the alpha particle orbits will be proportional to the 
angle between the alpha particle direction and the vertical axis, i.e., the neutron angle. 
Therefore, by measuring the alpha radius, one will measure the neutron direction for imaging. 
The radius can be measured in a number of ways known in the art. 

While the above detailed description has shown, described, and pointed out the 
fundamental novel features of the invention as applied to various embodiments, it will be 
understood that various omissions, substitutions, and changes in the form and details of the 
devices or processes illustrated may be made by those skilled in the art without departing from 
the spirit or essential characteristics of the invention. The described embodiments are to be 
considered in all respects only illustrative and not restrictive. The scope of the invention is, 
therefore, indicated by the appended claims rather than the foregoing description. All changes 
that come within the meaning and range of equivalence of the claims are to be embraced within 
their scope. 
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APPENDIX A 



A.l Gamma Energy 

The gamma energy peaks produced during the inelastic collision between fast neutrons 
and various atomic nuclei have a typical energy range of 0.100 to 18 MeV. The general 
contention in explosive detection has always been to detect the major components of explosives; 
that is, carbon, nitrogen, and oxygen. Therefore, the following discussion will focus on the 
gamma ray spectra of these constituent elements. It should be noted, however, that this principal 
is potentially applicable to many other elements including chlorine, phosphorus, and sulfur. 

Size and efficiency limitations of HPGDs limit the upper range of detectable gamma 
energies. For example, based on laboratory measurements with at 45% efficient N-type HPGD, 
a ceiling of 5.5 - 6.2 MeV prohibited the use of any gamma peaks about 6.129 MeV. The 
gamma energy peak selection process was therefore limited to gamma counts emanating below 
6.2 MeV. The 5.104 photo-peak of nitrogen produced 49 counts ± 14 using this detector. In 
comparison, an 81% N-type detector produced an energy ceiling of 7 MeV, with the equivalent 
5.104 nitrogen photo-peak (using equivalent neutrons) producing 197 counts ±21. 

While detector efficiency and size limits the ceiling of measured gamma peaks, the 
absorption of gamma rays through various materials influences the floor. Gamma attenuation is 
proportional to the thickness (Z) of interrogated materials. Further, based on measurements from 
experimental data of soil measurements, gamma rays below 1.6 MeV cannot be used for this 
application. Hence, the floor for gamma peak selection is at least in part based on the identity of 
interrogated material, and is on the order of 1.6 ± .5 MeV for most applications. 

A, 2 Cross Section 

After the determination of the gamma energy range, one must decide which gamma rays 
within the 1.6 to 7.2 MeV are to be used for analysis. One of the main factors for this selection 
process is atomic cross section. 

The probability of inducing a specific gamma ray depends on the binding forces within 
the nucleus of each atom. This probability is the atomic cross section or size (effective) for the 
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production of gamma rays and is measured in barns (10" cm). Established cross sections for 
carbon, nitrogen, and oxygen are used to establish a ratio between the number of carbon, 
nitrogen, and oxygen atoms. 

Carbon is somewhat unique with respect to the present analysis. First, carbon has a 
relatively simple gamma decay structure « there is only one detectable photo-peak (4.440 MeV) 
from 1.6 to 7.2 MeV. Second, the cross section of the 4.440 MeV peak (14 MeV incident 
neutrons) is 210.6 mb (very large). Third, the cross section of the 4.440 MeV peak increases 
dramatically as neutrons thermalize - 8 MeV incident neutrons have a cross section of 445 mb. 
The simple nature of the carbon decay scheme therefore translates the 4.440 MeV photo-peak to 
the ubiquitous carbon marker. However, the simple spectra of carbon has a consequence; 
because of the nature of the carbon nucleus, there is a recoil during the inelastic collision which 
produces a short-lived energy state. Subsequently, the photo-peak is rather wide with a shift of ± 
50 KeV. Thus, in order to obtain an accurate determination of carbon atoms, detection from 
4.340 to 4.540 MeV is needed. ' ■ . 

Nitrogen and oxygen have many photo-peaks from 1.6 to 7.2 MeV. The cross-section of 
the oxygen photo-peak at 6.129 MeV is 101 mb. Nitrogen has three detectable photo-peaks at 
2.312, 3.948, and 5.104 MeV with cross sections of 7.8, 3.5, and 31 mb, respectively. The 
minimum detectable cross section of a photo-peak is 5 mb (cascades often increase the 
observable cross section of photo-peaks, see section A.3 below). 

A3 Cascade vs. Photo-peak 

In addition to traditional energy levels, cascade effects increase the nurhber of possible 
energy peaks. Cascade effects are excitations of energy level that are not sufficient enough to 
cause a direct drop to the lowest energy state; and as a result, cascade peaks are produced from 
the systematic drop of photons from excited energy states to more stable ones. For example, the 
2.748 MeV oxygen cascade is produced by the excitation of the 8. ,872 MeV photo-peak and its 
subsequent drop to the 6.129 MeV energy level. The photon then drops from the 6.129 MeV 
(second photo-peak) to the ground state. In effect, the excitation of higher energy levels produce 
more gamma energy peaks and hence a greater probability of detection. 
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Unfortunately, data for cascade cross sections are not readily available. Experimental 
results have yielded approximate cross sections based on normalization to photo-peaks with 
detector efficiency and gamma absorption corrections. For oxygen, a normalized spectra of H2O 
elucidated an approximate cascade cross section of 28 mb for the 2.748 MeV peak. 

Cascade peaks represent over 50% of the total spectrum for oxygen and nitrogen. In 
order to maximize count rate, it is prudent to utilize cascade peaks to avoid the problems 
addressed by gamma energies outside the 1.6 to 7.2 MeV range, avoid the use of low cross 
section peaks, and finally circumvent the problem of proximity and overlap. 

A.4 Proximity and Overlap 

One of the most important factors in atometric analysis is proximity and overlap of 
gamma ray peaks. Proximity of gamma ray peaks is the simpler of the two cases. One of the 
advantages of using HPGe detectors is its increased resolution (on the order of 0.1%) at 622 
KeV. For instance, a HPGe detector can discriminate between a 5.156 MeV aluminum signal 
from a 5.104 MeV nitrogen signal. Nal detectors have resolution of roughly 10% at 722 KeV 
and cannot discern between many peaks in the spectrum. Proximity is an indication of the 
increased resolving power of the HPGD and is treated as an important tool for discrimination. 

Overlap is the extreme case of proximity. An example of overlap is witnessed in the 
carbon 4.440 MeV peak (recall the width of the peak is 100 KeV). Therefore, any gamma peak 
in this range of the spectrum would contribute to error in the carbon signal. The worst case is 
overlap by elements other than nitrogen and oxygen. Note that Aluminum has a photo-peak with 
a cross section of 4.9 mb at 4.41 1 MeV. . ♦ 

It is important to use information from the entire spectrum to scan for possible overlap 
and eliminate it by subtraction. If the overlap is between the three constituent elements, an 
appropriate algorithm may readily be developed to subtract extraneous counts. In conclusion, 
proximity and overlap determine the most 'opportune' gamma signals for use in the analysis. 

A.5 Single/Double Escape Peaks 
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Another important factor in determining which gamma energies are used for analysis is 
the production of single and double escape peaks. Pair production is the result of 
electron/positron annihilation with the concomitant release of 511 KeV (0.5MeV) quanta in 
roughly opposite directions. This process occurs in the matrix of the crystal lattice and is a 
function of the detector size. 

The most prominent single escape peak is from the oxygen 6.129 MeV peak. One sharp 
peak is produced 0.51 1 MeV lower at 5.618 MeV. A second escape peak is also produced 1.022 
MeV downstream at 5.107 MeV. Coincidentally, the nitrogen photo-peak at 5.104-5.107 peak is 
greater in intensity than the 5.618 MeV peak; this is a clear indication that there is a signal from 
nitrogen as well as oxygen. 

To eliminate the second oxygen escape peak from the nitrogen peak, one must determine 
the effective cross section of the escape peaks (specific to the chosen detector only). Based on 
experimental results of Si0 2 trials by Applicant, the cross section of the first and double escape 
peaks are 70% and 31% of the 6.129 MeV peak. In other words, the net number of counts in the 
5.107 MeV double escape, contributed from oxygen at 6.129 MeV, is 31% of net counts in 6.129 
MeV oxygen photo-peak. 
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APPENDIX B 



B.l Elimination of Noise 

A significant obstacle for quantitative atometry is the elimination of extraneous noise. N 
may arise from a variety of sources including air surrounding the test specimen and apparatus; in 
one m 3 of air, there is 925 grams of nitrogen and substantial oxygen. 

Accordingly, a method of dealing with this "noise" is needed. First, one must determine 
the volume of air interrogated by tagged neutrons. To accomplish this, a hypothetical 
rectangular box is used to represent the volume of interrogated air; For more precise 
measurements, one may replace the rectangle by a cone, to account for the drift of neutrons from 
the source. But for this investigation, the difference in volume is insignificant next to the size of 
the sample and an approximation with a rectangular volume will suffice. 

To determine exactly how much nitrogen and oxygen contribute to the data, it is 
necessary to calculate the dimensions of the rectangle. In one exemplary trial, the coincidence 
gate was open for a duration of 6.8 ns or roughly 34.5 cm centered on the target (neutron velocity 
is 5 cm per nanosecond). Since the target was exactly 1 meter away, the length of rectangle is 
34.5 cm — +/- 17 cm from the target. The area of the rectangle is a function of the alpha detector 
size. The alpha detector for this series of experiments had a measured area of 44.15 cm . The 
area of tagged neutrons, at 1 meter is exactly 4415 cm . Therefore, the total volume of 
interrogated air is [4415 cm 2 x 34.5 cm = 152,000 cm 3 ). 

The average density of dry air at room temperature is assumed to be p=0.00121 g/cm . 
As a result, it is possible to determine the mass of air [0.00121 g/cm x 152,000 = 178 g ]. Since 
the mass of the urea test sample is also known (5 kg in the present example) it is also possible to 
calculate the ratio of "moles of sample" to "moles of noise" (molecular weight of air ~ 29/gmole 
and urea 60.1 g/mole); specifically 6.87%. 

Finally, to determine the number of nitrogen and oxygen counts as a result of the air, the 
following formula was used (based on nitrogen): 

Percent nitrogen = 

(moles nitrogen air )/(moles nitrogen air + moles nitrogen urca ) x 100% (Eqn. B.l) 
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Equation B.l yields a correction value of 5.36% for nitrogen and 1.44% for oxygen. In other 
words, the net peak counts have to be corrected to account for extraneous signals from nitrogen 
and oxygen in air. 

It should be noted that the identity of each gamma line, to this point, has not been of 
consequence. However, based on the discussion in Appendix A above, any gamma line that 
corresponds to an overlap with an escape peak must be treated with special care. The nitrogen 
5.104 peak is this exception. An additional correction factor is required to calculate the 
percentage of oxygen (double escapes) that overlap with this signal 
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APPENDIX C 



C.l Correction for Detector Efficiency 

The concept of relative detector efficiency (i.e., that efficiency relative to a particular 
chosen energy value) is introduced herein. One of the reasons to use relative efficiency in the 
present invention instead of absolute efficiency is that only a quantitative atometric ratio of 
elements is needed, and not atomic content from one species outright. For instance, one does not 
have to solve implicitly for n (atomic density). Instead, if the ratio of atomic density is used, one 
may set C:N:0 equal to the ratio of n } to n 2 where the detector efficiencies are relative to each 
other. Table C.I lists experimental numerical values of the relative efficiency of a 80% N-type 
HPGe. 



TABLE CI 



Detector 


80% 


Relative to 


Emergy (KeV) 


Efficiency 


2300 KeV 


1800 


6.91E + 01 


0.14 


2300 


6.02E + 01 


1.00 


4440 


4.15E + 01 


1.46 


5105 


3.84E+01 


1.57 


6129 


3.46E+01 


1.74 



It should be noted that relative efficiency in Table C.I is based on the 2,300 KeV value. For 
example, the ratio of relative intensity increases 45% from 2,300 KeV to 4,400 KeV and 9% 
from 4,400 KeV to 5,100 KeV. 
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APPENDIX D 



D.l Background Subtraction - General 

It is evident that in order to quantitatively discern atomic information from gamma peaks, 
a method and algorithm for background subtraction be devised. The Regions of Interest ("ROI") 
for background subtraction depend on the identity of the interrogated material. In the case of 
explosives and drugs, these regions must overlap with the specific gamma signature lines of 
carbon, nitrogen and oxygen. A detailed list of gamma lines, escape peaks and cascades has 
been determined by the Applicant herein. 

D.2 Method 

First, the background on the upstream (lower energy) side of the peak is calculated as the 
average of the first three channels of the ROI. The channel number for this background point is 
the middle channel of the three points. Background on the downstream (higher energy) channel 
side of the peak is calculated as the average of the last three channels of the ROI. The channel 
number for this background point is also the middle channel of these latter three points. The two 
background points on each side of the peak form the end points of the straight-line background. 
Hence, the background is given by the following: 



( 1 + 1 h \ 



B = 



h-l + l 



IX+X C , (Eqn.D.l) 



where 



B — the background area 

1 = the ROI low limit 

h - the ROI high limit 

Cj = the contents of channel i 

6 = the number of data channels used (3 on each side in the present embodiment) 

The gross area (counts) is the sum of all the channels contained within the ROI according to the 
following: 
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A g =Y C, (Eqn. D.2) 



where 

A g = the gross counts in the ROI 
1 = the ROI low limit 
h = the ROI high limit 
Ci = contents of channel I 

The adjusted gross area is the sum of all of the channels within the ROI but not used in the 
background according to the following: 

A ag = g C, (Eqn. D.3) 

i=i+3 

where 

A ag = the gross counts in the ROI 
1 - the ROI low limit 
h = the ROI high limit 
Ci = contents of channel I 



B(h-l-5) ,r ^ 

The error in the net adjusted area is the square root of the sum of the squares of the error in the 
adjusted gross area and the weighted error of the adjusted background. The background error is 
weighted by the ration of the adjusted peak width to the number of the channels used to calculate 
the adjusted background. 
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APPENDIX E 



The following describes the overall method of the atometric process of the present 
invention. 

First, gamma information is gathered using an average value of peak width, location, and 
energy (Appendix A). Next, the gamma lines are systematically chosen for background 
subtraction (Appendix D). Each gamma line is subjected to numerous fitting routines until the 
lowest error value is repeated. Finally, corrections are made for extraneous noise, detector 
efficiency, and gamma attenuation (Appendix C). 

In order to ascertain the empirical formula for any contraband substance, the 
aforementioned data for one element are related to the number of gamma counts of another — 
such as carbon to oxygen — from the following established formula: 

I y -_F{dn\cjMyr ejr . rabs .v abs (Eqn. E.1) 

Here, I y = the gamma intensity = net gamma counts for specific energy range as calculated by 
background subtraction with the ROI algorithm. I T is standardized to a set number of channels; 
carbon =100, nitrogen and oxygen=12. r C fr, Tabs and v a bs are detector efficiency, gamma and 
neutron attenuation factors, respectively, and a is the cross section. Here, Feff includes solid 
angle subtended by the detector at 0.25 m (according to the relevant IEEE standard). 

The number of neutrons was calculated using the ratio of alpha to neutron counting rate: 



^ a \30neutrons ., „ 
fr = _ — .sec.- 1 — =# neutrons (Eqn. b.2) 

sec a 



F=neutron flux (number of neutrons into 4tt). 
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(8Q)-[7i(D 2 /4)]/(4ti 2 )-D 2 /1 6L 2 ) 



(Eqn. E.3) 



(5Q) = Solid r of object and L is distance from neutron source of the object. The number of 
atoms exposed to neutrons is: 



ft atoms — n.V = 



( atoms\( Volume^ 
V cm 3 A 1 J 



(Eqn. E.4) 



In this case, n = atomic density; and V = volume. And the cross product of atomic density and 
volume is atoms (n). 

Since atometric information is required from different elements in the same trial, it is possible to 
cancel out common factors. F, (5Q), and V abs . Equation E.l may be rewritten to cancel out 
factors in the following form: 



/ =a,nT eff .y 



eff V att 



(Eqn. E.5) 



Finally, Equation E.5 can be solved for n. 



n 



atoms 



oT eff . 



1 



7 



att 



(Eqn. E.6) 



Determination of the C:N:0 ratio is based on Equation E.6. yatten is a reciprocal because its 
effects decrease gamma value below 2,300 KeV and increase values above 2,300 KeV. E.6 may 
be written as follows: 



n 



r v att . 



atoms 



(Eqn. E.7) 



cffi 
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In order to solve of n^the noise- , efficiency - , and attenuation-corrected values (plus cross 
sections for carbon, nitrogen) may be used. For example, 



n ni ,ro g en = (76.5 x 1 .57)/(32 x 1 .3 1) (Eqn. E.8) 

= 2.06 



similarly, 



ricarbon = (256 X 1.24)/(210.6 x 1.24) (Eqn. E.9) 

= 1.04 



finally, 



n 0 xygcn = (127 x 1.41)/(101 x 1.74) (Eqn. E.10) 

= 1.02 



To reduce the atomic data into the empirical chemical formula, the atomic ratios are divided by 
the lowest value (1 .02) for oxygen. The complete empirical formula for urea is: 

C (1.04/1.02) = 1 (Eqn. E.ll) 

N (2.06/1.02) = 2.02 (Eqn. E. 12) 

0(1.02/1.02) = 1 (Eqn. E.13) 

By using the error values for the peak fitting algorithm and propagation of error, the complete 
formula of urea with error is as follows: 



Ci(±o9)N2.o2(+o.2i)Oi(±o.48) (Eqn. E.14) 
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APPENDIX F 



As previously described, the Dalitz Plot (or "A-plot") is a two dimensional graphical 
representation of three dimensions of information; namely, the carbon, nitrogen and oxygen ratio 
of explosives and other contraband. See Figs. 15a and 15b above. For example, the A-plot is a 
normalized atomic density for nitrogen, N'=N/(C+N+0), obtained from the gamma detectors, 
and represented by the perpendicular distance-squared from the horizontal side of the equilateral 
triangle. Similarly, the normalized atomic densities of carbon, C\ and oxygen, O', are 
represented by perpendicular distances-squared from the other two sides of the triangle. This 
information may easily be translated into an A-plot by the following steps (using urea as an 
example): 



Experimentally determined empirical formula = C1N2.02O1 



Carbon 



= 1 



Nitrogen 
Oxygen 



= 2.02 



= 1 



Total 



-4.02 



Partial Fraction C 



C/(C+N+0) 



(Eqn. F.l) 



C* = .248 



N' = .502 



O' = 248 



Square of partial fraction: 



(Eqn. F.2) 



C ,(2) = 0.0618 
N' (2) = 0.252 
0 ,(2) = 0.0618 



2 

Sum of squares (h ): 
h 2 = 0.376 



(Eqn. F.3) 



Square of each element over square of total = Y x 



(Eqn. F.4) 
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Y c = C 2 /C 2 +N 2 +0 2 = 0. 1 64 
Y n = 0.671 
Y 0 = 0.164 

Normalized square length = L c: (Eqn. F.5) 

L c = C' (2) /h 2 = 0.164 
L n = 0.671 
L Q = 0.164 

X and Y values of urea [107] in Dalitz Plot: (Eqn. F.6) 

X = .577350*(Y n +(2Y c )) = 0.57735 
Y = L n = 0.671 

(Eqn. F.7) 
(Eqn. F.8) 
(Eqn. F.9) 



Therefore, the boundaries of Dahtz Plot are: 
(0,0) 

(.57735,1) 
(1.1547,0) 
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